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La presente invention concerne une methode et un dispositif de detection multi- 
utilisateur. Plus particulierement, la presente invention concerne une methode et un 
dispositif de detection multi-utilisateur a maximum de vraisemblance pour un systeme 
de telecommunication DS-CDMA {Direct Sequence Code Division Multiple Access). 
5 Dans un systeme de telecommunication mobile DS-CDMA, la separation des 

communications en provenance ou a destination des differents utilisateurs est realisee 
en multipliant chaque symbole complexe d'un utilisateur par une sequence 
d'etalement qui propre a ce dernier, encore appelee pour cette raison signature de 
l'utilisateur. La frequence d'etalement (chip rate) etant superieure a la frequence des 

10 symboles, le signal transmis par chaque utilisateur est distribue (ou etale) dans 
l'espace des frequences. Le rapport entre la bande occupee par le signal etale et la 
bande occupee par le signal d' information est appele facteur d'etalement. A la 
reception, la separation d'un utilisateur donne est obtenue grace a un filtrage adapte a 
la signature correspondante. Lorsque le canal de transmission presente une pluralite de 

15 trajets de propagation, la sortie du filtrage adapte comporte autant de pics de 
correlation. Chaque trajet du canal peut etre modelise par un coefficient multiplicatif 
complexe et un retard. Les signaux s' etant propages selon les differents trajets peuvent 
etre alignes et combines au moyen de coefficients complexes conjugues des 
coefficients de trajet, realisant ainsi un filtrage adapte au canal de transmission. Pour 

20 simplifier la terminologie nous engloberons dans l'expression generate « filtrage 
adapte a l'utilisateur it » a la fois l'operation de filtrage adapte a la signature de 
l'utilisateur k et celle de filtrage adapte au canal de transmission. 

Pour combattre 1' interference entre signaux a destination (liaison descendante) 
ou en provenance (liaison montante) des differents utilisateurs, il a ete propose des 

25 methodes de detection multi-utUisateur et notamment des methodes de detection 
iteratives telles que celles connues sous le nom de PIC {Parallel Interference 
Cancellation) et SIC {Serial Interference Cancellation). Elles sont fondees sur 
l'iteration d'un cycle d'&imination des interferences comportant l'estimation des 
symboles emis, revaluation des interferences et leur soustraction des signaux recus. 

30 Bien que performantes, ces methodes ne sont pas optimales dans la mesure ou elles ne 
fournissent pas une estimation au sens du maximum de vraisemblance des symboles 
transmis par les differents utilisateurs. 

Une methode de detection multi-utilisateur a maximum de vraisemblance 
inspiree de l'algorithme de Viterbi a ete proposee par S. Verdu dans un article intitule 
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« Minimum probability of error for asynchronous Gaussian multiple access channels » 
publie dans IEEE Transactions on Information Theory, pages 85-96, Janvier 1986 
mais sa complexity est prohibitive car elle varie de maniere exponentielle avec le 
nombre d'utilisateurs . 

5 Plus recemment U a ete propose par L. Brunei et al., dans un article intitule 

« Euclidian space lattice decoding for joint detection in CDMA system » publie dans 
Proceedings of ITW, page 129, Juin 1999, une methode de detection multi-utilisateur 
a maximum de vraisemblance utilisant une representation par reseau de points. Selon 
cette methode on determine un vecteur caracteristique du signal recu representant une 

10 statistique suffisante pour la detection au maximum de vraisemblance des symboles 
transmis par les differents utilisateurs. On montre sous certaines conditions que le 
vecteur caracteristique peut etre represents comme le point d'un reseau perturbe par 
un bruit. La detection consiste alors a rechercher le point du reseau le plus proche du 
point correspondant au vecteur recu. Cependant, la dimension du reseau a utiliser etant 

15 de 2.K ou K est le nombre d'utilisateurs le nombre de points reste encore tres eleve. 
Pour simplifier la detection, il a ete propose de limiter la recherche du plus proche 
voisin aux points du reseau appartenant a une sphere centree autour du point re9u. 
Nous exposerons ci-apres cette methode de detection simplifiee dite « methode de 
detection par spheres » : 

20 Nous placons dans le contexte d'un systeme de telecommunication mobile a 

acces multiple avec etalement de spectre par sequences directes (DS-CDMA) 
comprenant K utilisateurs communiquant de maniere synchrone avec une station de 
base. 

Soit </*(/) le symbole complexe emis par l'utilisateur k a l'instant i . Ce symbole 
25 appartient a la constellation de modulation A k utilisee par l'utilisateur k, que l'on 
appellera encore alphabet de symboles de l'utilisateur k. 

Chaque utilisateur k transmet un bloc de N symboles avec une amplitude du 
signal ak. Les symboles sont etales par une signature complexe s k (t) = sf(t) + j.s k (t) 

de duree egale a la periode symbole T : 
30 ^(/) = 0si/«[0,r[ 



Les K symboles complexes dk(i)=dk(i)+j.dl(Q transmis a l'instant i sont places dans 
un vecteur ligne de valeurs reelles d2(/) defini comme : 



Du fait de la synchronisation des utilisateurs, le signal module correspondant peut 
s'ecrire, en fonction du temps / : 

&-^o*A(/>(/-iT) (2) 

i=0k=] 

Nous supposons que le canal est un canal ideal a bruit Wane additif gaussien. Soit 
5 n=St+?? le signal re9u au temps / et 7, un bruit gaussien complexe de moyenne nulle 
et dont les composantes ont une variance N 0 . 

Soit les vecteurs lignes yOMM'V et y^H^)^^--^)^)) oil 
yk(i)=yk(i)+ j.yl(i) est la sortie complexe a l'instant i du filtre adapts a Futilisateur k : 

A"* -00 

yt(f)= $sk\t-iT)ndt 

-00 

=^atdt(i)R£k+nk(i) 

T T 

10 avec Rtk=Ut(t)skXt)&t=R&+j.Rlk pour k,£ = l,.. m ,K et m(i)=fosk\t-i.T)di 

of 0 

On notera R(/) la matrice d'autocorr61ation des sequences d'&alement. 

Si Ton decompose les elements complexes de (3) en leurs parties reelles et 

imaginaires, on obtient : 



20 



15 \ym+JM0\=^ ( 4 > 
Soient A 2 =Diag(tfi,ai,.. et R 2 la matrice de taille 2Kx2K telle que : 



R2= 



1® Rir- Mc RIk 



-RkxRh—RhcRhc 



(5) 



L'equation (4) peut alors se mettre sous forme matricielle 



y2(/)=d2(/)M2+n2(i) 




ou M2 est une matrice reelle de taille 2K x 2K definie par M2=A 2 R 2 et ou le vecteur 
de bruit m(/H«,«(iV(4.^(i>Jt(0) a pour matrice de covariance N 0 R 2 . 

Nous demontrerons ci-apres que y 2 (/), tel que donne par l'equation (6), peut etre 
repr6sent<§ comme un point d'un reseau A 2 de dimension 2.K, de matrice generatrice 
5 M2 corrompu par un bruit m . 

Nous appellerons reseau reel de points A de dimension k tout ensemble de 
vecteurs de R* verifiant : 
x=6i vi +&V2+. . .+hc\K ou &eZ,V/=l v . .,ic 
et ou {vi,V2 y ..,v«-} est une base sur R* . 

Un exemple de reseau de points de dimension 2 a ete repr6sent6 en Fig. 1 . 

Les points du reseau forment un sous-groupe abelien additif de R K , c'est 

10 d'ailleurs le plus petit sous-groupe de R* contenant les vecteurs {vi,V2 r .,Vir} et un Z- 
module de R K . Ces vecteurs de base forment les lignes de la matrice generatrice G 
du reseau. On peut done ecrire x = bG ou h=4fn r .J>c)^L K ■ (7) 

La region delimitee par les vecteurs de base est appelee parallelotope 
fondamental et son volume, note vol(A) ou det(A) , est nomme volume fondamental. 

15 Ce volume fondamental n'est autre que le module du produit vectoriel des k vecteurs 
de base et est done egal a |det(G)| ou det designe le determinant. S'il existe plusieurs 

choix possibles pour la matrice g6neratrice d'un meme reseau, il n'existe par contre 
qu'une unique valeur pour le volume fondamental. 

La region de Voronoi V ou cellule de Dirichlet d'un point x appartenant au 
20 reseau est I'ensemble des points de R r plus proches de x que de tout autre point du 
reseau. Le volume de cette region est egal au volume fondamental. 

Le rayon d'empilement p du reseau est le rayon de la plus grande sphere 

inscrite dans la region de Voronoi et le rayon de recouvrement celui de la plus petite 
sphere circonscrite a cette meme region. Le rayon d'empUement est done le rayon des 

25 spheres dont l'empilement constitue le reseau de points et le rayon de recouvrement 
est celui des spheres les plus petites qui, centrees sur les points du reseau, permettent 
de recouvrir tout l'espace R* . La densite du reseau est le rapport entre le volume de la 
sphere de rayon p et le volume fondamental. Enfin, le coefficient d'erreur (kissing 
number) t(a) du reseau est le nombre de spheres tangentes a une m8me sphere dans 

30 l'empilement ou, en d'autres termes, le nombre de voisins d'un point du reseau, situes a 
la distance minimale = 2p . 
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Considerons a nouveau liquation (6). Les nombres complexes </*(/}f /.</*(/) 
appartiennent a un alphabet fini (ou constellation) A de cardinal : 

Card{A)=f\Card{Ak) (8) 

k=l 

Supposons par exemple que les composantes <//*(/) et dl(i) soient des symboles de 
5 modulation PAM d'ordre M : 

^(/>{-A/+l ,-M+3 r . ,M-3,M-l} et (9) 
<#(i)e{-M+1 . ^M-3,A/-l} (10) 

Si Ton effectue la transformation : 

d$ (i>l(^(i>A/-l) et diiiyfydlifYM-l) soit encore de maniere vectorieUe : 
15 d' 2 (/>4(d 2 (/>fv M ) (ID 



ou v« =(M-1 ,M-l , -M-\ ) 

les composantes d£(i) et <fi(/) sont des Elements de Z et par suite d^i) est un 

20 vecteur de Z 2IC . 

De maniere generate, s'il existe une transformation affine transformant les 
composantes «//(/) et dl(i) en des elements de Z, on peut representer le vecteur d'z^') 

par un vecteur de Z . 

De maniere similaire, on effectue la transformation correspondante sur y4f\ 

25 c'est-&-dire : 

yiO"H(y2(/>-VA/M2) (12) 



Moyennant cette transformation, que nous supposerons implicite dans la suite, le 
30 vecteur d2(0M 2 (i) appartient alors a un reseau de points A 2 de dimension 2.K tel que 
defini par la relation (7) avec G=M2(i) . Le vecteur y 2 (Q peut done etre considers 
comme un point du reseau A2 corrompu par un bruit 112(1)- 



Si Ton suppose que les composantes du vecteur de bruit n 2 (/) sont des variables 
aleatoires independantes gaussiennes centrees, le probleme de la detection au sens du 
maximum de vraisemblance des symboles emis par les differents utilisateurs se 
ramene a la recherche du point z 2 du reseau A 2 tel que sa distance a y 2 (z) soit 
5 minimale. 

En realite, les composantes du vecteur de bruit n 2 (/) sont correlees et la matrice 
de co variance de n 2 (i) est N 0 R 2 . 

Ann de se ramener au cas decorrele il feut efFectuer prealablement au decodage 
une operation de blanchiment du bruit. 

La matrice R etant hermitienne, la matrice d'autocorrelation R 2 est symetrique 
definie positive et peut done faire l'objet d'une fectorisation de Cholesky : 

R 2 =W 2 W 2 r < 13 > 
ou W2 est une matrice triangulaire inferieure de taille 2Kx2K. 

On definit un vecteur d'observation blanchie : y2(/)=y2(/)wr (14) 

ainsi qu'un nouveau reseau de points Q 2 constitue des vecteurs de composantes 
10 (^(/Xx/OX . jcHifcUif)) avec X2(i)=xJi)W^ , ou x 2 (/) est un vecteur de composantes 

(x?(i)jc{(f\. . ^cR(i\xlc(i)) appartenant a A 2 . 

On peut focOement montrer, qu'apres blanchiment, la matrice de covariance du 

bruit filtre mQJWf 1 est egale ZNohic ou I 2 * est la matrice identite de dimension IK. 

La detection comprend done une premiere etape de blanchiment du vecteur 
15 d'observation suivie d'une etape de recherche du plus proche voisin au sein du reseau 

de points Q 2 . 

Pour reduire le nombre de points a tester, comme illustre en Fig. 1, on peut 
limiter la recherche a une sphere centree autour du point y 2 . En pratique, le choix du 
rayon de la sphere resulte d'un compromis : U ne doit pas etre trop grand pour ne pas 

20 conduire a un nombre de points trop eleve et suffisamment grand pour inclure au 
mo ins le plus proche voisin. 

La Fig. 2 represente schematiquement un dispositif de detection multi-utilisateur 
utilisant une methode de detection par spheres. Le signal recu west fihree par une 
batterie de filtres adaptes a chacun des utilisateurs, 210,,...,210 K - Les composantes 

25 reelles et imaginaires du vecteur d'observation y 2 (/) en sortie des filtres adaptes sont 
transmis a une unite de calcul matriciel eflFectuant l'operation de blanchiment spectral 
selon la relation (14). Les composantes reelles et imaginaires du vecteur blanchi y 2 (0 
sont ensuite transmises a une unite de detection par spheres recherchant le plus proche 
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voisin du point re^u au sein du reseau Q 2 de dimension 2.K. Les coordonnees du plus 
proche voisin donnent directement les composantes reelles et imaginaires des 
symboles estimes dk(i) pour les differents utilisateurs. 

On peut montrer que la methode de detection par spheres telle qu'exposee ci- 
5 dessus possede une complexity en O(A^) ce qui peut severer tres p6nalisant lorsque le 
nocQbre d'utilisateurs est eleve. 
— / Le but de la presente invention est de proposer sous certaines conditions une 
simplification de la methode de detection par spheres. 

A cette fin, Pinvention est dSfinie par une methode de detection d'une plurality 
10 de symboles (d k (0) transmis par ou pour une plurality K d'utilisateurs, chaque 
symbole appartenant a une constellation de modulation et faisant l'objet d'un 
etalement spectral par une sequence d'etalement, ladite methode comprenant une 
etape de filtrage adapte pour fournir un vecteur complexe (y(*),y(i)) caracteristique 
dudit signal re?u, ledit vecteur complexe etant decompose en un premier vecteur 
1 5 ( y*(0.y*(0 ) et un second vecteur ( y '(i),y '(0 ) et au moins les plus proches voisins des 
premier et second vecteurs etant recherches au sein d'un reseau de points (A,Q) 
genere par lesdites constellations de modulation, les symboles transmis etant estimes a 
partir des composantes desdits plus proches voisins. 

Avantageusement, les sequences d'etalement (s*(0) sont constitu^s de multiples 
20 reels ( sl{i) ) d'un meme coefficient complexe (a) . 

Selon un mode de realisation de T invention, la recherche est limitee a un 
premier ensemble de points du reseau appartenant a une premiere zone predetermine 
(Z R ) autour du premier vecteur et a un second ensemble de points du reseau 
appartenant a une seconde zone pr6determinee autour du second vecteur. 
25 De meme la recherche peut etre Iimit6e a un premier ensemble de points du 

reseau appartenant a une premiere zone predeterminee (Zr) autour de Torigine et a un 
second ensemble de points du reseau appartenant a une seconde zone predeterminee 
(Si) autour de Torigine. 

Lesdites premiere et seconde zones predetermines sont par exemple des 
30 spheres. 

Avantageusement, la recherche du plus proche voisin du premier vecteur est 
effectuee sur une pluralite de composantes de ce dernier, la recherche etant limitee 
pour chacune desdites composantes k un intervalle defini par une borne inferieure et 
une borne superieure, lesdites bornes etant choisies de mantere a ce que ledit intervalle 
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ne comprenne pas des points relatifs a des symboles ne pouvant appartenir a la 
constellation de modulation. 

De meme la recherche du plus proche voisin du second vecteur peut etre 
effectuee sur une pluralite de composantes de ce dernier, la recherche etant limitee 
5 pour chacune desdites composantes a un intervalle defini par une borne inferieure et 
une borne superieure, lesdites bomes &ant choisies de maniere a ce que ledit intervalle 
ne comprenne pas des points relatifs a des symboles ne pouvant appartenir a la 
constellation de modulation. 

Prealablement a la recherche du plus proche voisin, le premier vecteur (y*(0) 

10 est avantageusement soumis a un traitement matriciel visant a substantieUement 
decorreler les differentes composantes de bruit de ce dernier. 

De meme, prealablement a la recherche du plus proche voisin, le second vecteur 
( y '(/)) peut etre soumis a un traitement matriciel visant a substantieUement d6correler 
les differentes composantes de bruit de ce dernier. 

15 Selon une variante de realisation de V invention, ladite etape de recherche est 

etendue a celle d'un premier ensemble de points plus proches voisins dudit premier 
vecteur, dits premiers voisins, et a un second ensemble de points plus proches voisins 
dudit second vecteur, dits seconds voisins, et les symboles transmis sont estimes de 
maniere souple a partir des symboles generant lesdits premiers et second voisins et des 

20 distances separant lesdits premiers voisins dudit premier vecteur d'une part et lesdits 
seconds voisins dudit second vecteur d'autre part. 

Selon un mode de realisation particulier de l'invention, on determine a partir des 
symboles estimes les contributions de chaque utilisateur aux signaux obtenus par 
l'etape de filtrage adapte et, pour un utilisateur k donn£, on elimine en sortie de l'etape 

25 de filtrage adapts les contributions des autres utilisateurs correspondant aux symboles 
deja estimes. Alternativement, on determine les contributions de chaque utilisateur au 
signal recu a partir des symboles estimes et, pour un utilisateur k donne, on elimine en 
entree de l'etape de filtrage adapte les contributions des autres utilisateurs 
correspondant aux symboles deja estimes. 

30 Si les symboles desdits K utilisateurs sont transmis de maniere synchrone, ledit 

reseau de points sera preferablement de dimension K. 

Si les symboles desdits K utilisateurs sont transmis de maniere asynchrone et se 
propagent selon une pluralite de trajets, la dimension du reseau sera preferablement 
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egale au nombre de symboles des differents utilisateurs pouvant interferer et non 
encore estimes. 

La pr£sente invention est egalement dSfinie par un dispositif de detection d'une 
pluralite de symboles (d k (/)) transmis par ou pour une pluralite K d 'utilisateurs, 
5 chaque symbole appartenant a une constellation de modulation et faisant Fobjet d'un 
etalement spectral par une sequence d'etalement, le dispositif comprenant des moyens 
pour mettre en oeuvre la m£thode exposee ci-dessus. 

Ce dispositif pourra en particulier utilise dans un recepteur de systeme de 
telecommunication mobile DS-CDMA . 
10 Les caracteristiques de Finvention mentionnees ci-dessus, ainsi que d'autres, 

apparaitront plus clairement a la lecture de la description suivante faite en relation 
avec les figures jointes, parmi lesquelles : 

La Fig. 1 represente un reseau de points utile a la methode de detection 
employee dans le recepteur illustre en Fig. 2; 
15 La Fig. 2 represente de maniere schematique la structure d'un recepteur DS- 

CDMA multi-utilisateur utilisant une methode de detection par spheres; 

La Fig. 3 represente de maniere schematique la structure d'un dispositif de 
detection multi-utilisateur selon un premier mode de realisation de T invention. 

La Fig. 4 represente de maniere schematique la structure d'un dispositif de 
20 detection multi-utilisateur selon un second mode de realisation de Finvention. 

L'idee a la base de Finvention est d'abaisser la dimension du reseaujde points en 
utilisant des sequences d'etalement de type particulier pour les difierents utilisateurs. 

Considerons a nouveau un systeme de telecommunication DS-CDMA avec K 
utilisateurs synchrones. Si nous choisissons des sequences d'etalement sk(t) a valeurs 
25 r^elles, les termes imaginaires de la matrice R 2 et par suite de la matrice M2(0 sont 
nuls. Par suite, on peut modeliser le systeme par un r6seau de points r£el A de 
dimension K et de matrice generatrice M(/) : 

y*(i)=d*(/)M(0+n*(0 (! 5) 

30 y'(z>d'(zOM(/)+ii'(0 (16) 

ou y R (i),d R (i),n R (i) (resp. y'^d'C/Xn'O) ) sont les vecteurs constitues des parties 
reelles (resp. des parties imaginaires) des composantes de y(0» d (0» n (0 » 
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M(/>AR(/*)ou R(i) est la matrice constitute par les coefficients fl«=Js/(/).»(/)d/ et A 

0 

est le vecteur des amplitudes des K utilisateurs. 

Les vecteurs d'observation y*(i) et y'(/) appartiennent a R*. Apres transformation 
eventuelle selon une relation du type de celle de (12), les vecteurs y*(/")et 
5 y'(/)peuvent etre consideres comme des points d'un reseau A de matrice generatrice 

M(/*) corrompus par du bruit. 

On montre aisement que les vecteurs de bruit n*(i) et n' (i) ont tous deux pour 
matrice de covariance M>.R(/) . R(0 etant une matrice symetrique definie positive on 
peut la factoriser selon une decomposition de Cholesky : R = WW r ou W est une 
10 matrice reelle triangulaire inferieure de taille KxK. Afin de decorreler les composantes 
de bruit, les vecteurs d'observation reels y*(i)et y'(/)sont tout d'abord soumis a une 

operation de blanchiment : 

y^)=yW (17) 
is yW(0w 7 "' ( 18 > 

Dans un deuxieme temps, on recherche les plus proches voisins des vecteurs 
y*(i) et y'(/') appartenant au reseau de points Q constitue des vecteurs 

xft^Ow 7 " 1 ou x(/) appartient a A. On peut montrer facilement, qu'apres 
20 blanchiment, les matrices de covariance des bruits filtres n K (ffW T et n'(i)W T sont 
toutes deux egales a M>Iat ou Ik est la matrice identite de dimension K. 

On voit done que le fait d'utiliser des signatures reelles conduit a une recherche 
de deux plus proches voisins dans un raeme reseau de dimension K alors que dans le 
cas general, complexe, le decodage necessite une recherche dans un reseau de 
25 dimension 2K. En fait on peut facilement montrer que le n§sultat precedent s'etend a 
tout jeu de signatures portees un meme nombre complexe, e'est-a-dire telles que : 
»(/)=cr.s*(0ou o est un nombre complexe et .y*(r)est reeL Ce sera bien entendu en 

particulier le cas si toutes les signatures sont imaginaires. 

La recherche du plus proche voisin pour le vecteur y'(i) se faisant selon le 
30 meme principe que pour y R (i), nous n'exposerons que la premiere. Cette recherche 
consiste a determiner le point x minimisant la metrique : 



■0 



u 
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/^V^lKH^Hf (19) 

1=1 

ou y*=x+n* et x=(xi,...,jc*r) est un point appartenant au reseau Q. 

Alternativement, on notera que le vecteur y*(/) n'a pas besoin d'etre blanchi si 
Ton utilise une metrique bas6e sur la matrice de covariance : 

nJiy R /x)=(y R -x)R-\y R -x) T (20) 

On designera par la suite, pour des raisons de simplification, z le vecteur 
d'observation blanchi (y*(i)) ou non (y*(0) et \\\\ la m&rique intervenant dans la 

10 relation (19) ou (20). 

Les points du reseau Q sont constituSs des vecteurs x tels que x=bG oil G est la 
matrice g&ieratrice du reseau et b=(fe^..,Ar), les composantes 6, appartiennent k 

l'anneau des entiers Z . 

Le d&ecteur restreint avantageusement son calcul de metrique aux points qui 
15 sont situes a Tinterieur d'une zone de la constellation situee autour du point re9u, 
preferentiellement a Pint&ieur d'une sphere de rayon donn6 Vc centre sur le point 
re?u z. Seuls les points du reseau situes k une distance quadratique inffcrieure k C du 
point re?u sont done consid6r6s pour la minimisation de la metrique. 

En pratique, le decodeur effectue la minimisation suivante : 

20 

^MI=™?JH (21) 

Pour ce faire, le decodeur recherche le plus petit vecteur w dans l'ensemble 
translate z-£l . On peut exprimer les vecteurs z et w comme suit : 

25 

z=pG avec pKPU-^PO 

w=^G avec S=(£,...^) ( 22 ) 

II est important de noter que p et £ sont des vecteurs r6els. Comme w=z-x oil 
30 x appartient au reseau Q, on a la relation £ = p k - 6, pour i=l^ . avec 
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Le vecteurw est un point du reseau dont les coordonnees £ sont exprimees 
dans le repere translate centre sur le point recu z. Le vecteurw appartient a une sphere 
de rayon quadratique C centree en 0 si : 



Dans le nouveau systeme de coordonnees defini par la sphere de rayon 
quadratique C centree en y est done transformee en une ellipse centree a Torigine. La 
factorisation de Cholesky de la matrice de Gram r=GG r donne r=AA r , ou A est 
une matrice triangulaire inferieure d'elements Sij . 

II faut noter que si le vecteur y a ete blanchi, il n'y a pas Ueu d'effectuer cette 
factorisation car la matrice generatrice de CI vaut AW et est done deja triangulaire 
inferieure. Toutefois, dans le cas ou Ton n'a pas precede au blanchiment prealable, la 
decomposition de Cholesky est necessaire. Dans tous les cas, on peut ecrire : 



H| 2 =Qfe) = ^GG r ^^C 



(23) 




(24) 



Enposant 




on obtient : 




(25) 



En s'interessant tout d'abord a la plage de variations possibles de & , puis en 
rajoutant une a une les composantes, on obtient les K inegalites suivantes, qui 
definissent tous les points a l'interieur de 1'ellipse : 
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qKK^K<C 

qK-\jc-igK-\+qK,K-\£icY+qKK£K 1 <C 



(26) 



On peut montrer que les inegalites (26) imposent aux composantes entieres de b de 
satisfaire a : 



Pk 



<b K < 




+ Pk-\ +<1kjc-\€k 



live 

-I 



C-<Ikk$k 



K ( K 
f-i+l \ v=*+l 



+p, + 



+ Pk-\ + 1kjc-\Zk 



(27) 



ou 



fx] est le plus petit entier sup^rieur au reel jcet [*J est le plus grand entier 
inferieur au reel x . 

10 Le decodeur possede K compteurs internes, k savoir un compteur par dimension, 

chaque compteur comptant entre une borne inferieure et superieure comme indiquS en 
(27), 6tant entendu qu'a chaque compteur est assocte un couple de bornes particulier. 
En pratique, ces bomes peuvent etre mises k jour de mantere recursive. 

Avantageusement, on liste toutes les valeurs du vecteur b pour lesquelles le 

15 point correspondant du reseau x = bG est situe en de<?& de la distance quadratique C 
du point re9U. Les points du n§seau situes en dehors de la sphere consid&6e ne sont 
pas testes. Avantageusement, pour chaque composante l 9 ..JC les bornes superieure et 
inferieure de la recherche seront ajustees ne mani&re a ne pas contenir des points qui 
sont certainement en dehors de la constellation. Ainsi les compteurs ne perdent pas de 

20 temps a parcourir des points qui, de toutes felons, ne sont pas des solutions. Par 




exemple, si tous les utUisateurs emploient une meme constellation de modulation 
PAM de taille Af, les bornes de recherche ne pourront sortir de Pintervalle [0,M-1]. 

De plus, la recherche a I'interieur de la sphere peut etre considerablement 
acceleree en mettant a jour le rayon Vc avec la derniere norme euclidienne calculee 
5 ||w||. Finalement, on selectionne comme meilleur point x celui associe a la plus petite 

norme ||wfl. 

De maniere a etre sflr que le decodeur trouve au moins un point du reseau, on 
choisit avantageusement un rayon de recherche superieur au rayon de recouvrement 
du reseau. On peut par exemple le prendre egal a la borne superieure de Rogers : 

10 

^^{KlogK+KloglogK+SKy^M (28) 

ou Vk est le volume d'une sphere de rayon unite dans l'espace reel . 

La Fig. 3 illustre schematiquement un dispositif de detection multi-utilisateur 
selon un premier mode de realisation de 1' invention. Le signal recu est tout d'abord 

15 filtre par une batterie de filtres adaptes aux differents utilisateurs, 310i,..,310 K . Le 
vecteur d'observation en sortie des filtres adaptes est decompose en un vecteur 
d'observation reel y"(0=0^(0»-^(0)et un vecteur d'observation imaginaire 
y'(i)=(yi(i),..,yK(i)). Apres transformation eventuelle du type de (12) (non 
representee), les vecteurs y R (/)et y'(i) subissent un blanchiment spectral en 320 et 

20 321 afin de decorreler les echantUlons de bruit. Les vecteur blanchis y*(i) et y'(i) font 
ensuite l'objet d'une recherche de plus proche voisin comme decrit plus haut dans les 
detecteurs par spheres 330 et 331. Le point trouve par le detecteur 330 donne 
(moyennant eventuellement une transformation inverse a celle de (12)) les 
composantes reelles des symboles estimes pour les K utilisateurs. De meme, le 

25 detecteur 33 1 donne les composantes imaginaires de ces symboles estimes. 

Au lieu de fournir directement des symboles de la constellation, le recepteur peut 
etre adapte pour fournir les symboles sous forme de decisions souples. Dans ce cas, la 
recherche a I'interieur de la sphere de detection ne se limite plus au plus proche voisin 
mais est etendue a une pluralite des plus proches voisins du point relatif au signal 

30 recu. 

Plus precisement, soit Z R et T,\ les spheres centrees respectivement autour de 
y*(i) et y'(/). On associe a tout couple d mM (Q=( vm,v*t) de points voisins appartenant 
a (Er,2i) tel que les K composantes <C"(0 appartiennent aux constellations de 
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modulation des utilisateurs, une probabilite a posteriori p mm , probability que le 
vecteur d? M (i) defini par ce point ait ete emis, etant donne l'observation yQ). Soit 0 
1'ensemble de ces couples. On definit un symbole souple d'un utilisateur k comme le 
A/*-uplet {m,..., 7tMk) ou Mk est le cardinal de la constellation de modulation de 
5 Putilisateur k et ou nj est la probabilite que le symbole s, ait ete emis. On a : 

^=pisj/y)= ^{sj/d^p^ (29) 

Les probabilites a posteriori p m ' m ' pourront par exemple etre exprimees en 
fonction des distances X m et Xm> separant les vecteurs y R (i) et y'(/) de V5„et v m >. 

L'equation (2) presupposait que les signaux des differents utilisateurs 6taient 
10 synchrones. Lorsque cette hypothese n'est pas valide, un symbole etale d'un 
utilisateur it aim instant donne peut interferer avec deux symboles Stales successife 
d'un autre utilisateur k\ Si nous supposons que la dispersion des retards de 
transmission t* des differents utilisateurs est inferieure a une p^riode symbole F, le 
symbole d'un utilisateur k emis a l'instant /, d k (i) peut interferer avec les symboles 
15 dk (i-l) et dk (i+l) d'un utilisateur k\ On peut supposer sans perte de generality que 
0<r\<.<Tk<.<TK<T . Aprds que les symboles d k (i-l) pour tous les utilisateurs k=J..K 
sont detectes, on commence la detection des symboles relatifs & l'instant i en 
commen9ant par l'utilisateur le plus pr^coce (ici l'utilisateur 1) et finissant par 
l'utilisateur le plus tardif (ici l'utilisateur K). La detection de d k (i) depend de trois 
20 vecteurs : le vecteur dp^i(0,...//*-i(0^ aux 
symboles deja detectes, le vecteur drK<ii(/+l),...^*-i(/+l),rf*(0,t/*+i(0, -A(0)relattf 
aux symboles futurs et le vecteur d'observation y(i) . On peut montrer que le vecteur 
d'observation complexe peut s'ecrire sous la forme d'une contribution passee et d'une 
contribution future : 

25 

y=dpARp+dFARp+n (30) 
ou n est un vecteur de bruit de matrice de covariance NoRf, Rp et Rf sont 
respectivement la matrice de correlation des signatures avec les signatures pass£es et 
la matrice de correlation des signatures avec les signatures futures. 
30 Lorsque les signatures sont £ valeurs reelles (ou plus g^neralement sont des 

multiples reels d'un meme nombre complexe) l'equation (30) peut se decomposer 
sous la forme de deux Equations ne faisant intervenir que des vecteurs reels, de 
maniere analogue aux equations (15) et (16) : 
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y H =dp R ARp+dp*ARF+n* ( 31 ) 
y'=dp / ARp+dp , ARF+n / ( 32 ) 

Au lieu d'effectuer la detection par spheres sur les vecteurs d'observation y* et 
5 y 7 celle-ci operera sur les vecteurs debarrasses des interferences dues aux symboles 
passes deja estimes a savoir : y R -dp R ARp et y'-dVARp 

Alternativement, au lieu de proceder a l'elimination soustractive des 
interferences au niveau du vecteur d'observation, cette elimination peut etre envisagee 
de maniere equivalents en amont, sur les signaux etales Les symboles estimes sont 
10 reetales spectralement et les contributions des differents utilisateurs au signal recu 
sont soustraites une a une, a la volee (au fur et a mesure de leur reetalement), en 
entree des nitres adaptes. Les entrees des filtres adaptes recevront done les 
signaux epures des contributions passees a savoir n-h(t) ou : 

15 /*(/)= Y^^/X/-/ T-n)+£a*dv(i-\ )»(/-(/ -1 ) T-n) (33) 
p<* *>* 

Le premier terme de droite dans l'equation (33) represente la contribution des 
utilisateurs qui sont en avance sur l'utilisateur k et pour lesquels par consequent une 
estimation des symboles </*•(/) est deja disponible. Le second terme represente la 

20 contribution des utilisateurs qui sont en retard sur l'utilisateur k: la contribution due 
aux symboles de(i) ne peut bien entendu etre evaluee a I'instant ou le symbole dk(i) est 
desetale. En revanche, la contribution due aux symboles precedents dn(i -1) peut etre 
evaluee puisque une estimation des symboles d*(i-\) est deja disponible. 

II importe de noter que, contrairement au cas synchrone ou l'etape de detection 

25 par spheres estimait simultanement tous les symboles dk(i), dans le cas asynchrone, 
les symboles des differents utilisateurs sont estimes les uns apres les autres suivant 
leur ordre d'arrivee. 

La Fig. 4 illustre un second mode de realisation de 1' invention mettant en oeuvre 
une elimination soustractive des interferences pour des utilisateurs asynchrones Pour 
30 simplifier la representation, seule la branche de detection concemant un utilisateur k a 
ete Ulustree. La somme I^t) symbolise la somme des contributions des autres 
utilisateurs et est soustraite en 405* au signal n en entree du filtre adapte 410*. La 
detection par spheres pour l'utilisateur k demarre a I'instant ou le symbole </*(/) est 
desetale. II faut noter qu'a cet instant les symboles interferents non encore estimes des 
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autres utilisateurs n'auront fah Pobjet que d'un d£setalement partiel (sauf bien 
entendu si les symboles interferents sont synchrones). Le module 430* (resp. 431*) 
op6re sur les composantes en sortie de 420* (resp. 421*). Alternativement les sorties de 
tous les modules 430* (resp. 431*) pourront etre demultiplexees pour etre traitees par 
5 un module commun 440 travaillant K fois plus rapidement. A la difference de la Fig. 
3, ici seule la Jfcieme sortie des modules 430* et 431* est utilise pour fournir le 
symbole estime </*(/). Le symbole <?*(/) est ensuite r6etale et la contribution 
akdk(i)sk(t-iT-n) de Putilisateur k est evaluee puis soustraite en entree de chacun des 

filtres adaptes 4 1 0* - via /* (t) . 
10 Si les canaux de transmission des differents utilisateurs sont du type multi-trajet, 

le probleme de relimination des interferences est plus complexe puisqu'il faut 
considerer les interferences entre tous les trajets de tous les utilisateurs. Le signal re9u 
peut s'ecrire : 

i=0 *=l /»=! 

ou Pk est le nombre de trajets du canal de transmission de Putilisateur k 9 ry>*=r*+ft>*est 
le retard cumute du retard (z*) k remission de Putilisateur k et du retard de 
propagation le long du trajet p du canal k et c P k est le facteur multiplicatif complexe 
associe a ce trajet. On supposera k nouveau que 0<n<..^Tk^.^TK<T et qu'en outre la 
20 dispersion des trajets est inferieure k la periode symbole : O<0 P k<T . II en resulte que 
0<r*<27\ 

Les filtres adaptes 410* effectuent alors un filtrage adapts a la signature de 
Putilisateur k et au canal de transmission k par combinaison MRC {Maximum Ratio 
Combining) des signaux des difiterents trajets. Plus precisement, le filtre 310 k effectue 
25 P operation suivante : 

y*(0=£cl* \sk\t-iT-T P k)n6t (35) 

Du fait de la dispersion des retards, la dur£e k prendre en compte pour 
30 Interference entre utilisateurs est de 2T et le symbole d'un utilisateur pourra 
interfftrer avec deux symboles consecutife non encore estimes d'un autre utilisateur. 
La detection conjointe porte alors sur tous les symboles interferents non encore 
estimes, les symboles interferents estimes etant, quant a eux, exploites pour 
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l'&imination soustractive. La complexity du ctetecteur par spheres est plus elev&s que 
dans le cas asynchrone mono-trajet puisque la dimension du r£seau a prendre en 
compte est de 2AT-1 au lieu de K. 

Enfin, il feut encore noter que l'estimation des contributions interferentes dans le 
5 cas asynchrone mono-trajet ou multi-trajet peut etre r£alis6e soit h partir de symboles 
de constellation (decisions dures) soit k partir des symboles souples des diflSrents 
utilisateurs. 

Bien que certains modes de realisation de Pinvention aient 6t6 represents sous 
forme de modules fonctionnels, il est clair que le dispositif selon l'invention peut Stre 
10 realise sous forme d'un processeur programme pour ex6cuter les diflKrentes fonctions 
representees ou sous la forme d'une pluralite de processeurs d£dies aptes a mettre en 
oeuvre une ou plusieurs de ces fonctions. 




REVINDICATIONS 

1) Methode de detection d'une pluralite de symboles (d k (0) transmis par ou 
pour une pluralite K d'utilisateurs, chaque symbole appartenant a une constellation 
de modulation et faisant Pobjet d'un etalement spectral par une sequence 
d'etalement, ladite methode comprenant une etape de filtrage adapte (310i,..,310/:) 

5 pour fournir un vecteur complexe (y(O.y(0) caracteristique dudit signal re9u, 
caracterisee en ce que ledit vecteur complexe est decompose en un premier vecteur 
(y*(i),y*(0) et un second vecteur (y'(0,y'(0) et qu'au moins les plus proches 
voisins des premier et second vecteurs sont recherches (330,331) au sein d'un reseau 
de points (A,Q) g£n£re par lesdites constellations de modulation, les symboles 
10 transmis etant estim^s a partir des composantes desdits plus proches voisins. 

2) Methode de detection selon la revendication 1, caracterisee en ce que les 
sequences d'etalement (**(/)) sont constitues de multiples reels ($?(/)) d'un meme 

coefficient complexe (a) . 

15 

3) Methode de detection selon la revendication 1 ou 2, caracterisee en ce que la 
recherche est limit6e a un premier ensemble de points du reseau appartenant a une 
premiere zone predeterminee (Sr) autour du premier vecteur et k un second 
ensemble de points du reseau appartenant a une seconde zone pr6d6termin6e (Ei) 

20 autour du second vecteur. 

4) Methode de detection selon la revendication 1 ou 2, caracterisee en ce que la 
recherche est limitee a un premier ensemble de points du reseau appartenant a une 
premiere zone predeterminee (S R ) autour de Porigine et a un second ensemble de 

25 points du reseau appartenant k une seconde zone predeterminee (Ei) autour de 
Porigine. 

5) Methode de detection selon la revendication 3 ou 4, caracterisee en ce que 
lesdites premiere et seconde zones predeterminees sont des spheres. 

30 

6) Methode de detection selon Pune des revendications precedentes, 
caracterisee en ce que la recherche du plus proche voisin du premier vecteur est 



# 
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effectuee sur une pluralite de composantes de ce dernier, la recherche etant limitee 
pour chacune desdites composantes a un intervalle defini par une borne inferieure et 
une borne superieure, lesdites bornes etant choisies de maniere a ce que ledit 
intervalle ne comprenne pas des points relatifs a des symboles ne pouvant appartenir 
5 a la constellation de modulation. 

7) Methode de detection selon Tune des revendications precedentes, 
caracterisee en ce que la recherche du plus proche voisin du second vecteur est 
effectuee sur une pluralite de composantes de ce dernier, la recherche etant limitee 
10 pour chacune desdites composantes a un intervalle defini par une borne inferieure et 
une borne superieure, lesdites bornes etant choisies de maniere a ce que ledit 
intervalle ne comprenne pas des points relatifs a des symboles ne pouvant appartenir 
a la constellation de modulation. 

15 8) Methode de detection selon Tune des revendications precedentes, 

caracterisee en ce que, prSalablement a la recherche du plus proche voisin, le premier 
vecteur (y*(i)) est soumis a un traitement matriciel (320) visant a substantiellement 
decorreler les diffSrentes composantes de bruit de ce dernier. 

20 9) Methode de detection selon Tune des revendications precedentes, 

caracterisee en ce que, prealablement a la recherche du plus proche voisin, le second 
vecteur (y'(0) est soumis a un traitement matriciel (321) visant a substantiellement 
decorreler les differentes composantes de bruit de ce dernier. 

25 10) Methode de detection selon Tune des revendications precedentes, 

caracterisee en ce que ladite etape de recherche est etendue a celle d'un premier 
ensemble de points plus proches voisins dudit premier vecteur, dits premiers voisins, 
et a un second ensemble de points plus proches voisins dudit second vecteur, dits 
seconds voisins, et que les symboles transmis sont estimes de maniere souple a partir 

30 des symboles generant lesdits premiers et second voisins et des distances separant 
lesdits premiers voisins dudit premier vecteur d'une part et lesdits seconds voisins 
dudit second vecteur d'autre part. 
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11) Methode de detection selon Tune des revendications precedentes, 
ca^acfefisee en ce que Ton determine a partir des symboles estimes les contributions 
"de chaque utilisateur aux signaux obtenus par l'etape de filtrage adapte et que, pour 
un utilisateur k donne, Ton elimine en sortie de l'etape de filtrage adapte les 

5 contributions des autres utilisateurs correspondant aux symboles deja estimes. 

12) Methode de detection selon Tune des revendications 1 a 10, caracterisee en 
ce que Ton determine (340) les contributions de chaque utilisateur au signal re?u a 
partir des symboles estim6s et que, pour un utilisateur k donne, Ton elimine en entree 

10 de l'etape de filtrage adapte les contributions des autres utilisateurs correspondant 
aux symboles deja estimes. 

13) Methode de detection selon Tune des revendications 1 a 10, caracterisee en 
ce que, les symboles desdits K utilisateurs etant transmis de maniere synchrony ledit 

15 reseau de points est de dimension K. 

14) M6thode de detection de selon la revendication 11 ou 12, caracterisee en 
ce que, les symboles desdits K utilisateurs etant transmis de maniere asynchrone et 
se propageant selon une pluralite de trajets, la dimension du reseau est egale au 

20 nombre de symboles des differents utilisateurs pouvant interfiSrer et non encore 
estimes. 

15) Dispositif de detection d'une pluralite de symboles (d k (/)) transmis par ou 
pour une pluralite K d'utilisateurs, chaque symbole appartenant a une constellation 

25 de modulation et faisant Fobjet d'un etalement spectral par une sequence 
d'etalement, le dispositif comprenant des moyens pour mettre en oeuvre la methode 
revendiquee selon Tune des revendications precedentes. 



16) Recepteur pour systeme de telecommunication mobile DS-CDMA 
30 comprenant un dispositif de detection selon la revendication 15. 
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